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We calculate the spectral density of the Dirac operator over an ensemble of configurations composed of overlap¬ 
ping instantons and anti-instantons. We find evidence that the spectral density diverges in the limit A —> 0. This 
indicates the breaking of chiral symmetry and also provides evidence that quenched QCD may be pathological in 
nature. 


1. Introduction 

The reader is referred to for background to 
this field. Preliminary results for this work have 
been given elsewhere |^. 

The QCD Lagrangian for Nf quarks in the 
chiral limit possesses a large family SUL^Nf) ® 
SUfi{Nf) ® Ub(1) ® Ua( 1) of global symmetries. 
Phenomenologically we find that the non U(l) 
part of this symmetry is dynamically broken to 
the diagonal subgroup, SUL{Nf) 0 SUii{Nf) 
SUoiNf)- This chiral symmetry breaking (ySB) 
overcomes the degeneracy between parity part¬ 
ners (particles with opposite parities but other¬ 
wise identical quantum numbers) and gives rise 
to N'j — 1 Goldstone bosons (a role played by the 
3 pions for Nf = 2). The order parameter for 
ySB is given by: 


= lim i 

m— 



2mp{X, m) 
+ m? 


( 1 ) 


where p{X,m) = limy^oo (p(A))ni. We note 
that the spectral density p(A,m) of the Dirac op¬ 
erator explicitly depends upon the mass m via the 
fermion determinant. Conventionally the above 
integral is calculated in the appropriate limits to 
give 


{iplp) = iTTp{0) 


( 2 ) 


This however can break down if the spectral 
density is not smooth in m. To illustrate, p 


= A““,q; G (0,1) implies (V'b) = hence 

the order parameter diverges in the chiral limit. 
If however p = (to/A)“ then we have a perfectly 
good order parameter in spite of the spectral 
density diverging in the limit of small eigenval¬ 
ues. Our first conclusion is therefore that one 
should use ( 0 ) to calculate the chiral condensate 
at various masses and then extrapolate to the chi¬ 
ral limit. One can however rely on the Banks- 
Casher relationship (i for quenched calculations, 
for here there is no fermion determinant to con¬ 
sider. 


2. Calculation of spectral density 

We construct an explicit matrix representa¬ 
tion of the Dirac operator i0[A] for a given 
gauge field A consisting of n_ instantons (I) 
and anti-instantons (I) (each object in iso¬ 
lation has gauge field Af). The matrix is 
of course defined by (i]p[A\)\j) = Dkj\k) in 
terms of some basis {|i)}. If the basis is or¬ 
thonormal then Dkj = {k\il/>\j). The basis we 
choose to use is the basis of “would be zero- 
modes” {1^:+),..., lbr)> • ■ ■: \'Pn.)} where 
{i]/)[Af])\il)f) = 0. We know that: 

= 0 

= 0 

= Vtj = yixt,pt,xj,p~) (3) 
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The first two follow from and 

{ 75 ,*-^} = 0. The function V (which is known 
approximately for a single I-I pair) in the third 
equation depends upon the center and size of 
the two instantons. In general it should also de¬ 
pend upon their relative colour orientation but 
we choose to ignore this, this should be irrelevant 
as long as there are no non-trivial correlation ef¬ 
fects in colour space. In terms of this basis we 
can construct the (n+ -|-n_) x (n+ -|-n_) matrix 
D = which has block zeroes on the di¬ 

agonal and V, respectively off block diagonal. 

There are two main objections to the matrix D 
being a matrix representation of the Dirac oper¬ 
ator. The first is that the “basis” we have cho¬ 
sen does not span the Hilbert space of wavefunc- 
tions. This is not too important for our consid¬ 
erations as wish to study the spectral density at 
small eigenvalues, these are the eigenvalues which 
have split from zero due to instanton interactions. 
The second more fundamental objection is that 
the “basis” we have chosen is not orthonormal 


(^+|^+) = U{x+,pt,x+,p+) 

= u{x-,p-,x-,pj) 

= 0 (4) 


hence the matrix D is simply not a representa¬ 
tion. We ignore the colour orientation of U for 
the reason given in the consideration of V. The 
third equation equals zero again due to the 75 
structure. We canjiowever construct a new or¬ 
thonormal basis {\ip~^), It/’”)} using the standard 
Gram-Schmidt procedure. 


1 < j < * < 

IV’7) = Slij\^~) l<j<i<n- (5) 

We are finally in a position to construct a ma¬ 
trix representation of the Dirac operator: 

ip = D = (tplipl'tp) 


where Dy = 


yt 0 


( 6 ) 


This matrix representation shares many of the 
properties of the Dirac operator (as it must). 


Firstly the matrix is endowed with a rich struc¬ 
ture, the elements are not random. This must 
be so as vectors in a Hilbert space obey triangle 
inequalities (a fact which is crucial is attempting 
to orthogonalize the basis). Secondly, non-zero 
eigenvalues occur in pairs ±A hence the matrix 
satisfies the 75 symmetry. Thirdly the Atiyah- 
Singer theorem is obeyed, a gauge field configu¬ 
ration with net winding number Q = |n_ — n+| 
does indeed have (at least) Q exact zero eigen¬ 
values. Fourthly our representation captures the 
essence of the mechanism by which eigenvalues 
are split away from zero. To see this simply con¬ 
sider the case of an I-I pair and notice that if 
we choose the function V appropriately then we 
will get the correct eigenvalue splitting. All these 
properties give us hope that the matrix represen¬ 
tation indeed shares the essential qualities of the 
actual QCD Dirac operator for determining the 
spectrum at small eigenvalues. 

The only other thing to decide is which func¬ 
tions to use for U and V ? Though these func¬ 
tions are known classically, there is little reason 
to believe that the classical forms dominate in the 
QCD vacuum. Pragmatically we use a variety of 
zero mode wavefunctions in calculating U and V 
(hard spheres, Gaussian and classical zero mode) 
and look for qualitative features of the resultant 
spectral densities which are independent of the 
functions used. In this talk we simplify and use 
the same functional form for both U and V, we 
discuss the general case (amongst other things) 
in a subsequent publication Q, but it sufhces to 
say that the conclusions are unaffected. 


3. Results 

Figure (H) shows the spectral density per unit 
volume for Nf = 0 (quenched) for two differ¬ 
ent volumes. The instanton gas is fairly dense, 
more dense than in the instanton liquid model 
Q but less so than found in some recent inves¬ 
tigations of the vacuum Q. We have used the 
classical zero mode wavefunctions for the over¬ 
lap integrals. The most noticeable thing about 
the resultant spectral density is the divergence 
at small eigenvalues. This divergence is indepen¬ 
dent of the functions used, see ||] for a more in 
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Figure 1. Classical Zero Mode V = U, (R^). Con¬ 
nected line from V = 97 -I- 97. Points (x) from 
V = 63 -I- 63, 100000 configurations. 


depth discussion. We can even use this method 
to generate a spectral density from real lattice 
data (where instanton positions and sizes are ex¬ 
tracted as a study of the vacuum). This is shown 
in figure (||). 

If we increase the density of instantons then 
we naively expect greater interactions, leading 
to eigenvalues splitting away further from zero. 
This should decrease the spectral density at small 
eigenvalues. This is indeed seen, as is the con¬ 
verse. The fact that we still see a divergence from 
the lattice data which is a dense gas shows the di¬ 
vergence to be a generic feature. Intriguingly we 
also see a difference if we position the objects at 
random, the spectral density increases at small 
eigenvalues. This shows that there may be non¬ 
trivial screening between instantons in the vac¬ 
uum. When we include the effect of the fermion 
determinant and simulate unquenched QCD then 
we still find a divergence. This however results 
in a perfectly acceptable quark condensate and 
{Q‘^)(m) distribution in the chiral limit. Our 
work however indicates that quenched QCD may 





Figure 2. Hard Sphere, (T^). (3 = 6.0, 32^64, 51 
configurations, (o) Positions and sizes from lat¬ 
tice data, (x) positions random. Data from Q]. 


be pathological in nature, it may have a diver¬ 
gent quark condensate. If this is so, then why 
has such a feature not been seen directly from 
the lattice ? We believe that the spectral den¬ 
sity is protected by lattice artefacts, there are no 
small eigenvalues on the lattice (unless by acci¬ 
dent), hence the spectral density cannot diverge 
as required. However such a feature should be 
apparent in the “ultimate quenched calculation”. 
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